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ABSTRACT 


An evaporative surface is described for heat pipes and 
other two-phase heat transfer applications that consists 
of a hybrid composition of V-grooves and capillary 
wicking. Characteristics of the surface include both 
a high heat transfer coefficient and high heat flux 


capability relative to conventional open-facet 
thread surfaces. With a groove density of 12. 


and ammonia working fluids heat transfer coefficients 
in the range of 1 to 2 W/cm^K have been measured along 
with maximum heat flux densities in excess of 20 W/crn^. 

A peak heat transfer coefficient in excess of 2.S W/cm^K 
was measured with a Z7.8 cm-1 hybrid stir face. 
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NOMENCLATURE 


a 

A, 

b 


c 




H 

H 



L 


M 


P 

p+ = 



Q 


a 


constant heat flux zone height, cm 

2 

cross-sectional area of single groove opening, cm 
bar height (Figure 2.3-1), cm 
friction factor coefficient (Table 2-1) 
local film thickness (Figure 5-2), cm 
hydraulic diameter of groove cross-section, cm 
laminar friction factor 

2 

gravitational acceleration, cm/sec 

fi-llet conductance, W/cmK 

latent heat of vaporization, ergs/gm 

2 

evaporative heat transfer coefficient, W/cm K 

2 

effective value of local film coefficient, W/cm K 

2 

local interfacial heat transfer coefficient, W/cm K 

wall and liquid thermal conductivity, respectively, W/cmK 

2 

capillary wick permeability, cm 
distance between vapor vent slots, cm 
working fluid molecular weight, gm/g-mole 
groove linear density, cm”^ 

covered groove perimeter, viewed in cross-section, cm 

ratio of differential vapor pressure to surface tension pressure 

2 

total pressure within vapor bubble, dynes/cm 

2 

liquid phase pressure, dynes/cm 

2 

saturated vapor pressure corresponding to T^^, .dynes/cm 

saturated vapor pressure corresponding to T_„, dynes/cm^ 

2 

vapor pressure within covered groove, dynes/cm 

2 

heat pipe vapor core pressure, dynes/cm 

2 

vapor pressure at liquid film surface, dynes/cm 

2 

evaporative heat flux density, W/cm 

2 

maximum evaporative heat flux density, W/cm 
heat flux per unit length of fi-llet, W/cm 

2 

local heat flux density at meniscus-vapor interface, W/cm 



V 


r = local film radius of curvature, cm 
r|^ = vapor bubble radius, cm 

rj^^. = vapor bubble departure radius from site with radius r^^. , cm 

r„, , = capillary pore radius for wick material covering grooves, cm 
cw 

r^^. = critical nucleation site radius such that ?■*■ = 1, cm 

r+. = r ./d 

m m 

R = radius of inscribed circle (Figure 2.1-1), cm 
Re = Reynolds number 

\h ~ thermal constriction resistance, cm°C/W 
t = capillary wick thickness, cm 

T = local wall temperature (Section 2.4), K - 
= mean liquid-vapor interfacial temperature, K 
= vapor phase temperature = K 
Tg = liquid film surface temperature, K 
= local wall temperature, K 
Tq = constant temperature at fin root, K 
T^r " liquid temperature at y = r, K 

Tsd = film surface temperature (saturated conditions assumed), K 
= liquid temperature at y = r, K 
= saturation temperature corresponding to K 
Twe “ exterior evaporator wall temperature, K 
= liquid film surface velocity, cm/sec 
= average vapor velocity in groove, cm/sec 
V^g = volumetric change per gram in vaporization, cmvgm 
w = fin half-width (Figures 2.3-1 and 2.4-1), cm 
w = land width (Figure 3-2), cm 

cl 

W|^ = groove width (Figure 3-2), cm 

We = Weber number 

X = coordinate defined in Figure 2.2-1 

= final height of vapor bubble if. bubble can grow from initial 
nucleation radius r^^.j, cm 



VI 


a 

e 

6 

Am 

AP 

AP 

AP 

AP 

n 

V 

e 

X 

u 


£in 


%o 


vm 


£vm 


£,v 


fi,,v 


cr 

<{> 


^ni = 


= defined by Equation 5-7b 

= fraction of local film thickness penetrated by vapor bubble, cm 
= local film thickness (equivalent to "d"), cm 

2 

= differential mass arrival rate at meniscus-vapor interface, gm/cm sec 

maximum liquid phase pressure drop relative to heat pipe vapor, 
dynes/cm^ 

liquid phase pressure drop up to position X = ^, relative to 
heat pipe vapor core, dynes/cm^ 

maximum vapor pressure rise in covered groove relative to heat 
pipe vapor core, dynes/cm2 

maximum liquid-vapor pressure difference, dynes/cm^ 

AP /AP„ 

vm £m 2 

liquid and vapor kinematic viscosities, cm^/sec 

total groove angle, radians 

characteristic dimension of surface disturbance, cm 
liquid and vapor phase viscosities, respectively, poise 

3 

liquid and vapor phase -densities, respectively, gm/cm 
surface tension, dynes/cm 

wall superheat factor defined by Equation 5-5, cm 
value of 4> corresponding to P"*" = 1, cm 

superheat factor corresponding to a bubble of height y = y^, cm 
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1.0 INTRODUCTION 

The development of methods for reducing the interfacial thermal resistance 
between a solid surface and a contacting liquid is an area of significant 
engineering and economic importance. One of the most effective techniques 
used is two-phase heat transfer which is characterized by a liquid-to-vapor 
transition at or near the solid surface. There are basically two methods of 
two-phase heat transfer, and this distinction relates to the predominant 
energy-carryi ng species. In boiling heat transfer, vapor bubbles generated 
at the surface create a turbulent mixing and pumping phenomenon that circu- 
lates unheated liquid from the bulk fluid into a position near the heated 
wall, and it is the sensible heat capacity in this colder liquid that, in 
most cases, dominates heat removal. In film evaporation, on the other hand, 
the surface is oriented or structured to minimize liquid film thicknesses, 
and heat removal is dominated by thermal conduction through the film and 
vaporization from the film surface. The heat pipe is a prime example of a 
device relying on film heat transfer. To minimize the film conduction tem- 
perature drop while maintaining uniform surface wetting, heat pipes have 
employed thin capillary wick surface coverings and grooves cut or extruded 
into the heat pipe wall. Over the past several years, capillary wick sur- 
faces have gradually lost favor because the grooved surfaces provide more 
reproducible performance and can be produced at low cost using extrusion or 
threading processes for groove densities less than about 50 cm"^. In addition, 

with groove densities on the order of 60 cm"^, heat transfer coefficients up 

2 f 11 

to 2 W/cm K have been measured with ammonia working fluid.' ‘ Pool boiling 

of ammonia, in- comparison, shows a heat transfer coefficient on the order of 

0.4 to 1.0 W/cm^K in the heat flux density range 1 to 10 For 

terrestrial applications, boiling heat transfer can be used in heat pipes, 

but for space applications, boiling in zero-g appears to produce a large vapor 

blanket around the heated object which severely limits heat transfer. 

Disadvantages of grooved -surfaces for film heat transfer include poor capil- 
lary pumping capability and fabrication difficulties when more than 50 grooves/ 
cm are desired. A grooved surface has a periodic liquid film structure with 
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a principal radius of curvature equal to approximately one-half the groove 
width, while a capillary wick has two principal radii of curvature. For the 
same characteristic liquid fillet dimension, a capillary wick will pump twice 
the static head or twice the heat flux (assuming identical axial heat pipe 
fluid return conductances). 

High groove densities are desired because it is commonly held that the majority 

of heat transfer occurs over a relatively short zone where the liquid meniscus 

attaches to the groove. When conventional machining processes are used 

to produce over 50 grooves/cm, difficulties arise in uniformity of groove 

shape and depth, chip removal, and tool wear. In the production of axial 

(7) 

grooves, a similar limit occurs in the range of-40 grooves/cm. ' ' 

The hybrid surface to be discussed here is a combination of capillary wick 
and grooves that possesses the pumping capability of the capillary wick along, 
with the high heat transfer coefficient of fillet evaporation. 
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2.0 ANALYSIS 

2.1 Description 

A cross-sectional view of the evaporative surface is presented in Figure 2.1- 
la, while Figure 2.1-lb shows a corresponding section through a V-groove sur- 
face. The hybrid surface, which has been coined the "inverted meniscus wick," 
is produced by covering the V-grooves with a fine-pored capillary screen or 
felt-metal. The- capill ary material does not cover the entire surface, however 
The coverins is periodically, vented perpendicular to the groove axes to allow 
the grooves to communicate with the general vapor space. Prior to heat addi- 
tion, the capillary wick and grooves will be primed with liquid, but when the 
metal substrate is heated, liquid will recede into the groove because of wick 
pressure drops and perhaps as a result of boiling caused by wall superheat. 

The resulting liquid/vapor structure will be similar to that shown in Figure 
2.1-la, with -an inverted liquid meniscus and a tear-drop shaped vapor micro- 
channel. Vapor generated through thin-fi1m evaporation from the fillet 
structure passes down the vapor microchannel and is exhausted into the general 
vapor space through one of the periodic vent slots.. In contrast to open 
groove surfaces, allowable fillet recession is now totally independent of 
groove size and shape, and instead, is tied to liquid pressure drops and 
capillary wick pore size. Application of this evaporative wick to a cylindri- 
cal spiral groove evaporator is .shown in .Figure 3-3. In this case, capillary 
screen wick was wrapped on a steel mandrel, and the mandrel and wick were then 
inserted -into the grooved tube. Vapor generated under the screen cover 
traveled a circumferential path in the array of spiral grooves, and was dis- 
charged into the general heat pipe vapor core through two axial slots in the 
mandrel . 

As depicted in this figure, the inverted meniscus wick is more complex- than 
the open groove evaporator, and perhaps more costly (although attention has 
not been specifically directed towards mechanical design optimization). 
Experimental tests described. in a later section show the inverted meniscus 
possesses very high .heat transfer coefficients and maximum heat flux densities 
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and these attributes .must be balanced against the increased surface complexity 
to determine whether the inverted meniscus surface is appropriate for a specific 
application. 

2.2 Operating Constraints 

Liquid and Vapor Pressure Drops 

The operating limits of an inverted meniscus evaporator deal primarily with 
hydrodynamic pressure drops, and particularly with ensuring that the vapor 
phase pressure drop in the groove does not exceed capillary wick pumping capa- 
bility. In a previous report, pressure drop constraints were derived for 
V-groove surfaces. In this section, that analysis will be broadened to in- 
clude grooves of arbitrary shape. 

Figure 2.2-1 illustrates two fundamental vapor-liquid flow geometries for the 
inverted meniscus. The vapor slots are shown spaced a distance "L" apart. 

In Figure 2.2-la, the flow of liquid and vapor are counter-current, while in 
Figure 2.2-lb, the liquid and vapor flows are co-current. The first condition 
corresponds to a situation where the vapor exhaust slots at X = ±L/2 coincide 
with the points of liquid introduction, while the second condition corresponds 
to a situation in which liquid is introduced at the mid-point, X = 0, For 
condition -la, It is apparent that the maximum vapor/liquid pressure difference 
occurs at X = 0, and is the algebraic sum of the liquid pressure change in 
the capillary wick and the vapor pressure change in a vapor microchannel . The 
liquid pressure contribution is 

^2 

8 ' ^pA^gt (dynes/cm^) (2.2-1) 

where is the liquid-vapor pressure difference at X = ±~, and where the 
common assumption of uniform evaporation per unit has been employed. 
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FIGURE 2.2-1. TWO ALTERNATE METHODS OF FEEDING LIQUID TO THE SURFACE WHICH RESULT IN COUNTERCURRENT AND COCURRENT 
LIQUID AND VAPOR FLOW AT THE EVAPORATOR WALL 
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Estimation of the vapor phase pressure difference between X = 0 and either 
X = ±L/2 requires calculation of the microchannel hydraulic diameter Dj^, and 
the cross-sectional area Ag. For laminar vapor flow and uniform evaporation 
per unit length of meniscus, the vapor pressure difference between the mid- 
point and the vapor release slot is given by 


AP 




(dynes/cm^) (2,2-2) 


Use of the following formula has been made: 

P = -£- 

Re 


(2.2-3) 


where F^^ is the laminar friction factor, and the coefficient "c" is given in 
Table 2-1. For a triangular groove with total groove angle 0, the hydraulic 
diameter is 


D = • ■ 9' = 
^h P 


W[j cos (0/2) 

' 1 + sin (0/2)J 


(2.2-4) 


The total liquid-vapor pressure difference for configuration 2.2-la is the 
sum of equations 2.2-1 and 2.2-2. For the special case that Ap^^ = 0, a 
maximum heat flux density q^„ is defined on the basis of available capillary 

ain 

pumping pressure^ 



^am^ ^v 

8p£^gt 


1 , o/^v\ ct 

K ^ '■\v I 2 

V \z/ NAgD^ 


(2.2-5) 


The peak heat flux density q,_ and corresponding vapor flow Reynolds number 

aiH 

are shown in Figure 2.2-2 as a function of ammonia saturated vapor temperature 
for the experimental 37.8/cm grooved surface described in Section 3. In the 
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TABLE 2-1. LAMINAR FRICTION 
(F = c/Re) 

FACTOR coefficient! 

Covered Groove Cross-Section 

Coefficient "c" 

tube (ref.) 

16.0 

rectangle (3.5:1) 

17.75 

square 

14.25 

equilateral triangle 

13.25 

isosceles right triangle 

13. ,00 


1 


From H. Schlicting, Boundary Layer Theory , McGraw-Hill, New York, 1968. 


range -10*^ to -20°C, the maximum heat flux density is calculated to be 20 to 
30 W/cm^, while at 25^0, the maximum heat flux density is about 45 W/zwr. 

This increase is primarily caused by a lower vapor phase pressure drop, which 
in turn is the result of an increased vapor density. The conclusion to be 
drawn is that quite high heat fluxes are possible before vapor core pressure 
drop is significant. Because of the large number of small parallel vapor 
channels, vapor flow Reynolds numbers on the order of 500 or less are typical 
for covered V-grooves with high-pressure fluids like ammonia. 

For the case in which the liquid and vapor flow are co-current (Figure 2.2-lb), 
the maximum vapor-liquid pressure difference occurs generally at an inter- 
mediate point 0 < X < L/2 because the highest vapor pressure is found at X = 0, 
while the lowest liquid pressure is found at X = ±L/2. Therefore, the largest 
difference between these two pressures in general occurs at a position between 
X = 0 and X = ±L/2. 
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If the liquid loss rate by evaporation from the capillary wick is assumed to 
be uniform over the area of contact with the wall, and if the rate of vapor 
injection into each microchannel is assumed independent of position, then the 
liquid and vapor pressures are, respectively, 



and 


P 


V 



The largest difference between (2.2-6) and (2.2-7) occurs at 


(2.2-7) 


X _ 1 

L/2 ■ 1 + n 


where 

q 






( 2 . 2 - 8 ) 


and the maximum pressure difference is 




n 

1 + q' 


(2.2-9) 


Examination of this equation shows that co-current flow 'reduces the maximum 
vapor-liquid pressure difference. For example, for the particular situation 
that q = 1 (APyj^ = AP^) and AP^^ = 0, the maximum vapor-liquid pressure 
difference at the wall will be. 
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(countercurrent flow. Fig. 2a) 

(2,2-10) 


(cocurrent flow, Fig. 2b) 

(2.2-11) 


The peak heat flux density would therefore be about 33% higher for co- 
current vapor-liquid flow. 

Cocurrent flow may also be of benefit from a quite different standpoint. If 
liquid entrainment occurred in the vapor stream, this entrained liquid would 
be transported toward the groove end that was in most need of liquid, and 
nearest to a dry-out condition, whereas entrainment in countercurrent flow 
would accelerate dry-out. The subject of entrainment is discussed in more 
detail in the following section. 

Entrainment 

A second phenomenon that may occur with wick-covered grooves is liquid entrain- 
ment in the groove vapor stream or liquid wave formation because of liquid/ 
vapor interfacial 'shear. These effects are both correlated by the Weber 
number. 


Pv ^ 

We = ^ (2.2-12) 

CT 

where X is the characteristic wavelength of the disturbance. 

To a degree, entrainment and wave formation may be desirable because of poten- 
tial improvement in heat transfer resulting from turbulent .wave motion and 

from entrained droplet impingement on non-wetted groove surfaces. For a 

f 91 

Weber number = 2 tt, large interfacial waves and entrainment are expected.' ' 

The characteristic wavelength A for covered grooves will be on the order of 
the inscribed diameter 2- R (Figure 2.1-1). Weber numbers based on operating 
data are presented in Figure 2.2-3 for groove densities of. 12. 6 cm~^ and 37.8”^ 
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FIGURE 2.2-3. MICRO-CHANNEL WEBER NUMBER FOR EN- 
TRAINMENT BASED ON INSCRIBED CIRCLE OF 
RADIUS "R" 
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(with groove characteristic cross sections as in Figure 3-2), and it can be 
concluded that heat flux densities in excess of 50 W/cm are required for the 
onset of entrainment or wave formation. In actual fact, the critical heat 
fluxes are even higher because the operating menisci are significantly smaller 
than 2 • R at high power levels. 

2.3 The Role of Thermal Constriction 

Feldman and Berger^^^ and Schneider et al.^^^^ have found that two-dimensional 
finite element conduction models of meniscus heat transfer indicate that most 
evaporative heat transfer from a grooved surface occurs over a relatively 
small zone near the fillet tip where the liquid film is quite thin. The heat 
transfer coefficient for a grooved surface is therefore dominated by heat 
transfer over a very small fraction of the total groove area. The overall 
thermal resistance of this small area can be considered the sum of a geometri- 
cal constriction resistance in the underlying metal plus a film conduction 
drop across the meniscus plus an interfacial resistance at the liquid-vapor 
interface. Of these three resistances, the liquid-vapor interfacial resistance 
is generally small, and will be neglected at this. point. However, the rela- 
tive importance of constriction and film conduction in defining evaporative 
heat transfer behavior is not as obvious, and a simple model of constriction 
heat transfer for rectangular grooves was generated. 

In Figure 2.3-1, heat is removed at a uniform rate per unit area from an 

u 

infinitely long bar over a strip of height "a". The heat is input to the bar 
over the bottom face, which is at a constant temperature. All other areas 
are insulated. The detailed calculations will not be presented, but the 
general solution technique was as follows. The heat conduction equation was 
solved for the bar using separation of variables and the following boundary 
conditions. 


B.C.l - at X = 0, T = Tq 

B.C.2 - at y = w, 3T/9y- = 0 

B.C.3 - at X = b, 8T/3X = 0 
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■> 


INSULATED 

BOUNDARY 



FIGURE 2.3-1. CONSTRICTION RESISTANCE IN AN INFINITELY LONG BAR OF HEIGHT b AND 
WIDTH w. ONE SIDE IS AT A CONSTANT TEMPERATURE, WHILE HEAT IS 
BEING APPLIED ON AN ADJACENT SIDE OVER A STRIP OF HEIGHT "a", AS 
SHOWN. 
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B.C.4 -aty = 0, b-a, 9T/3y = 0 

B.C.5 - at y = 0, b-a rs x ^ b, 3T/3y = -Q /K 

3 W 

where T 1s the base temperature and is a constant heat flux density. 

Once the temperature distribution was found, the average temperature over the 
spot b-a < X ^ b was calculated. This immediately yielded the constriction 
resistance as 




O c» 

” l✓_-l 


K=1 


coth 




(K-%)' 


. cos^ [it{K-%){l-|)] (2.3-1) 


/Ol- 

where is the bar thermal conductivity and is the heat flux per unit 
length of bar (watts/cm). 

The constriction factor is presented in Figure 2.3-2 as a function of 

a/b for various ratios of w/b. The constriction resistance for a constant 
temperature, spot boundary condition was found for a specific value of w/b 
using graphite resistance paper. A comparison of these resistance values 
with the constant heat flux model (equation 2.3-1) is given in Table 2.3-1. 
For small spot, sizes, the constant temperature constriction resistance is 
about 6 % higher than the constant heat flux. However, it must be recognized 
that the accuracy of values derived from graphite resistance paper measure- 
ments is uncertain, and this percent difference cannot be totally ascribed 
to the differing boundary conditions. 


For constriction resistance to be of significance in limiting film heat trans- 
fer, it must be shown to be comparable in value to the evaporative resistances 
demonstrated for various grooved heat pipes. A great deal of data exists for 
axial square-grooved and circumferential V-grooved heat pipes with groove 
densities on the order of 25/cm to 35/cm. For example, Schlitt, Kirkpatrick, 
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TABLE 2.3-1- A COMPARISON' OF THEORETICAL CONSTRICTION RESISTANCE VALUES 
FOR THE INFINITE STRIP ^FIGURE 2-3-T) WITH VALUES OBTAINED 
■ USING GRAPHITE RESISTANCE PAPERH ) 


a/b 

Thermal 

Resistance 

(Analog) 

Thermal 
Resistance 
(Fourier series) 

Rf 

4.3100"^) 

7.65 

7.09 

1.08 

6.08(10“^) 

7.31 

6.87 

1.06 

— i 

CO 

o 

1 

CO 

6.97 

6.446 

1 .08 

1.5300"^) 

6.72 

6.28 

1.07 

7.3100"^) 

5.50 

5.29 

1.04 

1.2250 0"^) 

5.12 

4.96 

1,03 

.0344 

4.53 

4.30 

1.05 

.0569 

4.19 

3.977 

1.05 

.0732 

4.03 

3.813 

1.06 

.122 

3.-635 

3.473 

1.05 


^ ^The analytical model incorporates a constant heat flux at one boundary, 
whereas resistance paper necessarily approximates constant temperature 
boundary conditions, w/b = 0.3235 
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FIGURE 2.3-2. CONSTRICTION RESISTANCE FOR HEAT TRANSFER TO A RECTANGULAR BAR VIA A NARROW STRIP OF 
WIDTH "a" 
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and Brennan measured a heat transfer coefficient of 0.7 W/cm^K with a 8.7 
axial groove/cm aluminum heat pipe using ammonia at If the 

heat transfer coefficient of Schlitt, et al. is considered, to be the result 
of a series of parallel line heat transfer elements (that is, the meniscus 
tips), then the conductance per unit element length, or per unit fillet 
length, is 


h " 2 x*8 7 " 0-0402 W/cmK 

From Figure 2.3-2, the constriction resistance of a square groove (w/b = 

0.50) is about 5.7 at a/fa = 10 and approximately 7.2 at a/b = 10"^. Tak- 
ing the most pessimistic case of a/b = 10~^, this corresponds, for a 8.7/cm 
square groove, to a constriction width a = 5.75 x 10 cm, and a per-unit 
length fillet conductance of 

% " ~ ” 0.328 W/cmK (aluminum) 

This hypothetical contact zone of 575 Angstroms width has a lineal conductance 
that is a factor of 8.1 higher than the measured value. However, for stain- 
less steel, which has a conductivity on the order of 0.15 W/cmK, the lineal 
conductance would be only 0.0208 W/cmK, and constriction resistance may be 
of some importance. For high conductivity metals, the wall-side conduction 
drop does not appear to be a dominant factor in defining the overall evapora- 
tive heat transfer coefficient. ' 
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2.4 One-d1mensionaT Meniscus Heat Transfer Model 

To gain a fundamental understanding of meniscus heat transfer, a one-dimensional 
evaporative heat transfer model has been developed that not only accounts for 
conduction in the fillet and underlying wall, but also considers fluid dynamics 
within the fillet. Two-dimensional models of fillet heat transfer generated 
by other investigators^^~^^ have indicated that heat transfer primarily occurs 
over a very small region near the fillet tip. This means that a significant 
liquid flow must occur between each fillet base end tip, and that this flow 
must be pumped by surface tension forces. Since surface tension forces are 
intimately related to liquid-vapor interfacial curvature, it is quite possible 
that significant changes in surface curvature occur between the base and tip 
of each meniscus, leading' to a flow-influenced film curvature and heat trans- 
fer coefficient. 


For high thermal conductivity metals like aluminum, and for rectangular grooves, 
the constriction resistance analysis of Section 2.3 indicates that the con- 
duction drop through the heat pipe wall and ground surface tooth is signifi- 
cantly smaller than the conduction drop across the liquid film, on the basis 
of comparisons with experimental evaporative heat transfer data. Therefore, 
a one-dimensional "fin equation" approach to heat transfer in the evaporator 
wall is considered adequate for engineering model purposes. 


Fillet Hydrodynamics 


If the liquid fillet is viewed as a hydrodynamic boundary layer, the governing 
Navier-Stqkes equation is 


3U. 


3x 


+ U. 


3U, 


3y 


= v„ 


2 

9 U 


9y 


j_!!i 

3x 


(2.4-1) 


where the coordinate system originates from the fillet tip and x- is parallel 
to the groove wall and y- is perpendicular to the wall (Figure 2.4-1). The 
liquid kinematic viscosity is and the liquid density is p^. The liquid 
pressure P, is assumed independent of y. 
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In addition to the x-directed Navier-Stokes equation, conservation of mass 
gives 



3l = 

9y 


(2.4-2). 


The coupled system of equations (2.4-1) and (2.4-2) has as unknowns, U^, U^, 
and 1^. To reduce this system to one more amenable to solution, the von 
Karman-Pohl hausen integral technique has been used. In essence, an educated 
guess is made of the. form of the velocity distribution U . The boundary layer 

A 

equation (2.4-1) and conservation equation (2.4-2) are integrated at constant 
X across the liquid film. The result is a more tractable coupled set of 
ordinary differential equations that satisfies (2.4-1) and (2.4-2) in an 
average sense. For the present problem the x-directed velocity is assumed to 
be parabolic in y. 


= A + By + Cy^ (2.4-3) 

At the wall, the x-directed fluid velocity is zero, while at the liquid-vapor 
interface the effects of shear are neglected (3U /3y = 0). The coefficients 
in equation (2.4-3), with these boundary conditions, yield the velocity profile 

' ^1 - <a)^) <2.4-4) 

where 5 is the local film thickness and is the x-directed liquid velocity 
at the liquid-vapor interface. The y-difected velocity can now be obtained 
from the continuity equation as 

% - -f -w < 2 -'‘- 5 > 

0 

where use has been made of the boundary condition = 0 at y = 0. 
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When the expressions (2.4-4) and (2.4-5) defining and respectively , are 
substituted into the Navier-Stokes equation, and the equation is Integrated 
over the range 0 :s y < 6, the following equation in and 6 results: 



„2 ^ 
dx ■ 15 ‘^l dx 


S dx 


(2.4-6) 


To obtain a second equation, it is necessary to consider the relationship 
between film thickness and evaporation or condensation mass transfer. Specif- 
ically, for a small zone as shown in Figure 2.4-2, mass conservation requires 
that 



(2.4-7) 


where is the local heat flux density and h^g is the latent heat of vapori- 
zation per unit mass. The heat flux Q, is positive for evaporation and nega- 

a 

tive for condensation. In the limiting case of small Ax and A5, and using 
equation (2.4-5), 



3 ^a 


(2.4-8) 


Equations (2.4-6) and (2.4-8) would now be solvable if Q_ and dP/dx were known. 

a 



§8 
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§ 





FIGURE 2.4-1. SCHEMATIC DRAWING OF MENISCUS HEAT 
TRANSFER MODEL SHOWING COORDINATE 
SYSTEM AND CHARACTERISTIC DIMENSIONS 
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FIGURE 2.4-2. MASS BALANCE AT LIQUID-VAPOR INTERFACE. THE 
SYMBOL "d" IS EQUIVALENT TO THE SYMBOL "fi" IN 
THE TEXT. 
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Heat Flux Model 

The local heat flux is assumed to be determined by a film heat transfer 
coefficient (K^/S) and a vaporization limit. The rate at which molecules 
undergo a phase change at the liquid-vapor interface is finite and governed 
by kinetic theory. Under the assumption that each molecule impacting on 
the surface is thermalized to the film temperature (accommodation coefficient 
= 1.0), the differential arrival rate is 



(gm/cm^ sec) (2.4-9) 


where M is the liquid molecular weight and and respectively the 

film surface and vapor phase pressures. The corresponding film surface and 
vapor phase temperatures are T^ and T^^- The factor R is the gas constant. 

If the working fluid saturation vapor pressure is given by an equation of 
the form 

Pvs " ^ (2.4-10) 


where P^, N and are empirical curve-fit factors, and if the film vaporiza- 
tion temperature drop is small compared to the average system temperature T^, 
then the interfacial heat transfer coefficient is given by 







(N + h) 


) 


(W/cm^K) (2.4-11) 


This vaporization heat transfer coefficient and the film coefficient form 
a series thermal network with an overall heat transfer coefficient H from 
wall to vapor phase. 
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(2,4-12) 


where 



Xn 

£ m 


■33 hfg Vm - (N,J5) tJ 


(cm) (2.4-13) 


and where is in atmospheres, h^^ is in Joules/g-mole, is in K, and 
the liquid conductivity is in W/cm-K. The working fluid molecular weight 
is in gm/gm-mole. For ammonia at 20°C, the film vaporization process is 

thermally equivalent to the conduction drop across an ammonia film of 

-6 

1.5 X lO” cm thickness. 

The heat flux density Q_ is then given by 

a 

Qa = IT (T^j- T^,). Wa/) (2.4-14) 

where T^^ and T^ are the vapor phase and wall temperatures. 

Pressure Gradient Model 

For momentum transfer to be described by a unidirectional boundary layer 
approximation to the Navier-Stokes equations, the pressure gradient ,9P/8y 
in the liquid film must be assumed everywhere to be negligibly small. The 
only significant pressure change occurs in the direction of flow; that is., 
in the x-direction. If the fillet has only a one-dimensional curvature, 

-and the vapor phase above the fillet is assumed to be at constant pressure, 
then changes in liquid pressure within the fillet must correspond exactly 
to changes in surface curvature: that is, the changing surface curvature 

is actually the mechanism by which liquid is fed to the meniscus tip. This 
situation is shown in Figure 2.4-1, The liquid pressure at any film cross- 
section is virtually constant, and hence the pressure difference at vertical 
position "x" between the point "£1" in the liquid and point ”vl" in the vapor 
is 


ORIGINAL PAGB H 
OP POOR QTJAL^ 
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(2.4-15a) 


where r^ is the local one-dimensional liquid surface curvature and a is the 
liquid surface tension. At a position slightly further down the film, 

"X + Ax", the corresponding pressure balance is 



{2.4-15b) 


Since the vapor pressure is assumed constant limit of Ax 

becoming vanishingly small the liquid pressure P depends on the curvature "r" 
in the following way. 

4 1 (?) 

Using the cartesian coordinate expression for surface curvature, the liquid 
pressure drop is 



Conduction Drop in Metal Tooth or Fin 

For modeling of heat transfer through the underlying metal tooth or fin, the 
classic one-dimensional fin equation has been used. The metal tooth of 
Figure 2 is assumed to be a fin of width 2w, and the local fin temperature 
Ty is given by 


d^T. 


dx 


w _ H 


wk ^^w " ^vs^ 
w 


(2.4-16) 


where is the metal conductivity, and H is given by Equation 2.4-12. 
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Model Summary 

By incorporating the expressions for 
final coupled differential equations 
velocity and wall temperature 



5 * 6 ” 15 6 * dx 


Q^,dP„/dx, and T, developed above, the 

Q. W 

in film thickness 6, the film surface 
are 



(2.4-17) 


dUi o Q • /i\ U- .j, 

= . 1 (1\ J, . ^ (2.4-18) 

dx 2 \6/ 6 dx 


and 


d^T. 


dx^ 


w _ 


wk 


w 


- Tvs) 


(2.4-19) 


The resulting equations are third-order, first-order, and second order. The 
equations 'were solved using a coupled 4th order Runge-Kutta integration method. 
If the equations are posed as follows. 


y’" = f(y, y‘, y", u, T) (t = Tj 


u’ = g(U, y, y' , T) 

ORIGINAL B 

, , OF POOR QUMiW 

T" = h(T, y) 


(2.4-20) 

(2.4-21) 

(2.4-22) 


then the integration formulas are as follows, where "i+l" denotes the 
Runge-Kutta predicted value at the next step and "i" represents the present 
value of the variable. 



2-25 


^ ^ ^ ^ '^4^ (2.4-23) 

^i+1 " ^ ^ *^4^ (2.4-24) 

Ay2 

y;+l = y; + ^ (4 + k^ + k 3 ) + Axy-' (2.4-25) 

^i+1 " ^i ^ ^ ^ *^2^ ^ (2.4-26) 

■^i+1 = 2 N 2 + -2N3 +-N^) . (2.4-27) 

Ax2 

"^i+1 " ^i ^ T "^2 ■*■ ^3^ (2.4-28) 


To obtain the first values K^j W^, and N^, the equations 20 through 22 are 
evaluated at the immediate position "i"; that is 

= f(y^-, yi, ri, U., T.) (2.4-29) 

Wj = g(u^., y., yl, T.) - (2.4-30) 

= h(T., y.) (2.4-31) 

To obtain the values K 2 , W 2 , and N 2 > the equations 20 through 22 are evaluated 
with the following variable values. 


y = ’^i'^T^i (2.4-32) 

Ay 

y' = yi+x^i' (2.4-33) 

y" = ■*■ T “^1 (2.4-34) 

u = u^. + ^ (2.4-35) 

T = T. + f Tl (2.4-36) 
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To obtain the values W^, and the equations 2,4-20 through 2.4-22 are 
evaluated with the following variable values. 

2 

y = y^- +^y: (2.4-37] 

2 

y' = y;-+^y-+^K^ • (2.4-38] 

y" = y- + X *^2 (2.4-39] 

U. = u. +“W 2 (2.4-40] 


T - T + T' + A.^ . M 

I *i ^ 2 'i 4 


To obtain the values K^, W^, and N^, the equations 2.4-20 through 2.4-22 are 
evaluated with the following variable values. 

2 3 

y = y^ + Axy^l ^ ^ '^1 (2.4-42) 

Ax^ 

y' = y' + AxyV ’+ K 2 (2.4-43) 


y" = yV + AxK^ 


U. + AxW^ 


T. + AxT; + ^ N 2 


Although the problem statement appears formidable and the third-order 
differential equation has a variety of non-linear terms, the solutions were 
found to be quite well-behaved. 
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Starting the Solution 

Runge-Kutta integration is basically a shooting technique, and for a coupled 
third order differential equation system, initial values must be assigned for 
6, d5/dx, d 5/dx , U^, T, and dT/dx. For reasons of computing economy, it is 
preferable to begin the solution at the fillet tip, x = 0. At x = 0, the 
intuitive substitution 6 = 0 is not satisfactory since several terms in the 
differential equations are divided by 6. Instead, 6 must be set to a value 
less than the equivalent film evaporation thickness, which in the case of 
ammonia is about 1.5 x lO"*^ cm at 20°C. This initial value does not then 
dominate heat flow at the meniscus tip, and a local "sonic limit" is therefore 
assumed. The initial value of d6/dx is the tangent of the contact angle, 

A contact angle of 1° has been assumed for initial calculations. The second 
derivative d 6/dx is undetermined and can be arbitrarily set. However, it 
has been found that physically realistic meniscus conditions dictate a rather 
narrow range of appropriate second derivatives. If, for- example, the initial 
value of d 5/dx is too small, the film will not grow properly away from the 
meniscus tip, and in some cases the film will actually thin down and vanish, 
much like the simulation of an evaporating droplet on a hot surface. The 
value of T (the wall temperature) at the tip is also undetermined, but the 
insulated behavior for x = 0 requires dT/dx = 0 at x = 0. At the boundary 
X = 0, the x-directed film surface velocity has been set to zero. In summary, 
typical starting conditions for ammonia modeling are; 

1. x = 0 

2. 6=1 (10"^) cm 

3. d6/dx = 0.01746 

4. d^6/dx^ = variable 

5. = 0 

6. T = variable 

7. dT/dx = 0 

The extremely low equivalent film conduction drop for evaporating ammonia 
(1.5 X 10"° cm equivalent film thickness) poses a physical interpretation 
problem. For an actual film thickness much less than this value (1.0 x lO"^ 
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cm, for example), the molecular structure of the fluid may render this 
continuum analysis invalid. At a film thickness of 1.0 x cm, the film 
is only a few molecules thick. 

This problem also shows up in calculated meniscus profiles as a need for an 

2 2 

extremely high starting value for d 6/dx . The heat flux density is so high' 
in this region that the predicted curvature pressure is on the order of an 
atmosphere below the general vapor core pressure. This means, among other 
things, that there could actually be a film boiling process at the meniscus 
tip, not attributable in total to wall superheat, but also to a very low 
liquid phase pressure. 

Sample Solution 

Figures 2.4-3a and 2.4-3b present results computed for an ammonia meniscus 
in contact with a rectangular aluminum tooth of 0.020 cm {the width "w" is 
therefore 0.010 cm). In Figure 2,4-3a, the film thickness and percent of 
total heat transfer are shown as a function of distance from the meniscus tip. 
Total meniscus depth is about 0.0105 cm, measured from the meniscus tip to 
the meniscus centerline position,, while about 50% of the total heat is trans- 
ferred in the first 10% of this distance, confirming the importance of this 
thin-film region. 

In Figure 2.4-3b, the film thickness calculated with the present analysis is 

compared with the film thickness that would be expected if the meniscus were 

cylindrical in shape, as has been generally assumed for conduction modeling 

-3 

of meniscus heat transfer. For distances less than 0.5 (10 ) cm, this 

comparison has not been done in Figure 2.4-3b because of the non-zero contact 
angle. However, for distances greater than 0.5 (10 ) cm, where contact 

angle effects are negligible, the predicted film thickness is from 1.5 to 5 
times greater than the cylindrical meniscus approximation. This is caused 
by a liquid phase pressure drop parallel to the underlying wall that is 
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FIGURE 2.4-3 a. h. TYPICAL COMPUTATIONAL RESULTS FOR A SINGLE MENISCUS 
CONTACTING A METAL TOOTH OF 0.020 CM TOTAL WIDTH. 

IN FIGURE 2.4-3b, THE COMPUTED FILM THICKNESS IS COMPARED 
WITH THE THICKNESS EXPECTED FOR A CYLINDRICAL MENISCUS 
OF RADIUS Fjj. BECAUSE OF LIQUID FILM HYDRODYNAMICS, THE 
FILM THICKNESS DEVIATES CONSIDERABLY FROM THAT 
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associated with liquid transport to the meniscus tip.- This analysis indi- 
cates, at least for the limiting case of meniscus recession to the wetting 
angle, that fillet hydrodynamics is a crucial parameter in evaporation from 
grooved surfaces. 

If a meniscus conductance is defined on the basis of the temperature difference 

between the wall at x=x^ and the saturated vapor, then the lineal fillet 

conductance is 0.068 W/cmK. If it is assumed that the actual physical system 

corresponds to a 0.020 cm square groove, then the conductance based on the 

groove root temperature is 0.066 W/cmK. For a 0.01°C superheat at the fillet 

tip, the lineal heat flux is about 0.0173 watt/cm, which translates into 
2 

0.865 watt/cm for a uniform groove density of 25/cm, and a heat transfer 
coefficient of 3.3 W/cm^K. 

Conventional V-grooved ammonia/ aluminum heat pipes with 25 grooves/cm have 

2 

typically shown heat transfer coefficients, on the order of 0.5 to 0,8 W/cm K, 

while data presented in a later section for the inverted meniscus show a heat 

transfer coefficient of about comparable magnitude for a 37.8/cm grooved surface 
2 

at 0.865 W/cra . Therefore, the heat transfer coefficient predicted by the 
present model is about a factor of 5 higher than actual measured values. Addi- 
tional analysis is required to determine the basis for this rather large dif- 
ference. On a speculative level, this difference may be attributable to larger 
contact angles, nucleate boiling within the meniscus, or to a two-dimensional 
constriction resistance in the underlying wall, although the latter is not con- 
sidered highly probable, as discussed in Section 2.3. 

In other calculations for groove densities on the order of 10/cm, fillet con- 
ductances of about 0.075 W/cmK were calculated, while the axial U-groove data 
of Schlitt et al.^^^^ for 8.7 grooves/cm, as discussed in Section 2.3, yield a 
fillet conductance of 0.0402 W/cmK. Experimental fillet conductances-, there- 
fore, appear to be in better agreement with calculations for U-groove surfaces 
with low linear groove densities. 
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3,0 EXPERIMENTAL EQUIPMENT 

To test the performance of the inverted meniscus concept with various groove 
densities, a heat pipe was constructed with a short cylindrical evaporator 
section as detailed in Figure 3.1. The active evaporator surface is 3.81 cm . 
long and 1.27 cm in diameter for a total internal surface area of 15.20 cm . 
The aluminum evaporator assemblies were fabricated with circumferential 
groove densities of 12.6 cm”^, 25 cm”^, and 37.8 cm"^, as in Figure 3.2. The 
grooves in each evaporator were cut as a single continuous spiral groove with 
a nominal 65° included angle. After the grooving process, a finishing cut was 
taken along the axis of the evaporator to produce the characteristic flat- 
topped groove shape shown in Figure 3.2. This final cut produces a uniform 
internal tube diameter and minimizes the possibility of groove tip deformation 
during installation of the liquid transport wick. The 25 cm"^ evaporator was 
constructed for an earlier study, and different machining techniques re- 
sulted in a groove profile resembling a half-cylinder rather than a V-groove. 

The evaporator heater in each case consisted of a single layer of spiral - 
wrapped Inconel ribbon. The ribbon material was 0.16 cm wide, and was wrapped 
on a single layer of 0.009-cm-thick Teflon tape. The gap between each spiral 
wrap of ribbon was 0.11 cm, and was used as a pass-through zone for evaporator 
thermocouples. A layer of DC 340 applied to the heater winding assisted in 
maintaining uniform and adequate thermal contact of the ribbon with the Teflon 
film. The ribbon heater assembly and thermocouple leads were potted in a 
flexible silicone blanket such that the blanket has an outside diameter 
identical to the outside diameter of the evaporator end flanges. The silicone 
material stabilized the relative position of the heater ribbon and thermo- 
couple leads so that the possibility for direct contact of the heater ribbon 
with the thermocouple leads was minimized. 

The evaporators were each thermally monitored with 12 thermocouples inserted 
radially to a depth of 0.46 cm in the 0.62 cm thick evaporator wall. Eight of 
the thermocouples were inserted at 90° intervals in proximity to the two end 
flanges, while the remaining four thermocouples were spaced at 90° intervals 
in the center of the evaporator. 
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FIGURE 3-1. 




V-GROOVE THREADS 



X 

1.266 

1.269 


EVAPORATOR ASSEMBLY FOR INVERTED MENISCUS TESTING IN CYLINDRICAL 
GEOMETRY. ALL DIMENSIONS IN CENTIMETERS. 
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CHARACTERISTIC VALUES 

DIMENSION (CM) GROOVE DENSITY (GROOVES/CM) 




12.6 

25.0 

37.8 

Wa 


0.0260 

0.0081 

0.0069 

Wb 


0.0534 

0.0319 

0.0196 

hg 


0.0421 

0.0161 

0.0151 

R 


0.0146 

0.0078 

0.0053 


FIGURE 3-2. CROSS-SECTIONAL VIEWS OF SPIRAL GROOVE 
EVAPORATING SURFACES 


7705 - 124.10 


3-4 


An adiabatic section 45 cm in length was used for all tests, while a con- 
denser length of 15 cm was used for the 25 cm~ tests, and a 30 cm length for 
the tests performed with 12.6 and 37.8 cm~^ grooved evaporators. A blind- 
end, thin-walled stainless steel tube 0.32 cm in diameter and 30 cm in length 
was inserted from the far end of the condenser to monitor vapor core tempera- 
ture. 

A close-fitting guide tube covered with seven layers of 400 mesh stainless steel 

screen served as a wick to distribute liquid circumferentially to the screw 

-1 -1 
threads for the 12.6 cm tests, while six layers were used in the 25 cm 

tests. A cross-sectional detail of this wicking structure is presented in 
Figure 3.3. The guide tube incorporated two 0.108-cm-wide axial slots so 
that the screens may be drawn tightly over the tube and inserted through the 
slots to be spot-welded to the inner surface of the guide tube, and to pro- 
vide a means for vapor escape from the grooves. The slots together account 
for 5.4% of the interface area. The layers of screen are allowed to overlap 
on the inside guide tube diameter to provide an area for the axial wick to ' 
transfer fluid to the circumferential wick. As assembled, the screen-covered 
tube is perhaps .0.0025 cm larger than the I.D. of the grooved evaporator, but 
inherent looseness in a multilayer wick allows the guide tube and wick to be 
inserted. 

To ensure adequate fluid transport to the evaporator section, the 'axial wick 
shown in Figure 3.4 was used. This multiple-tube artery consisted of 19 
screened single-layer tubes of 0.076 cm O.D. fitted into an outer wrap of 
0.457 cm diameter. All screens used were 400 mesh square-weave stainless 
steel. The liquid was transported from the axial wick to the circumferential 
wick by four multiple-layer screen tubes. In the 25 cm"^ tests, the tubes 
consisted of 6 layers of screen, while later tests with 12.6 and 37.8 cm~^ 
surfaces used 8 layers of screen. Except for very low power operation, the 
central holes in these screened tubes remained unprimed. A summary of heat 
pipe mechanical design and zero-g transport capacities are given in Table 3.1. 
Axial flow conductance of the artery was calculated using the hydraulic 
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FIGURE 3-4. SECTIONAL VIEW OF INVERTED MENISCUS HEAT PIPE EVAPORATOR WITH 

MULTIPLE-TUBE ARTERY AND SCREENED TUBE RADIAL TRANSPORT WICKS. 
ALL DIMENSIONS IN CM. 
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Parameter 


TABLE 3.1. HEAT PIPE SPECIFICATIONS 
12.6 and 37.8 cm“^ Tests 


25.cm“^ Tests 


Evaporator I.D.^^^, cm 
Evaporator length, cm 
Adiabatic length, cm 
Condenser length, cm 
Axial transport wick 


Max. capillary pumping « 
pressure (exptl.), dynes/cm'^ 


Calculated transport 
capacity, watts @ 25°C^^^ 


1.27 
3. BI- 
AS. 1 
30.5 


1-.27 

3.81 

45.1 

15.2 


19 tube bundled artery, 
400 mesh square weave screen 


12,600 at = 20.1 dynes/cm 
772/407 545 


TT5 ! 

^To inner faces of threaded section 
( 2 ) 

' ^Includes vapor flow pressure drop 
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diameter method discussed in Reference 11, while liquid and vapor pressure 
drops in the evaporator were calculated following the procedures discussed 
in Section 2.1. Capacity predictions corrected for Ig but not for a liquid 
pool, are shown in Figure 3.5, along with various test points. Test point 
adverse tilts are corrected for liquid puddling, but the overfill used to 
ensure priming of the multiple tube artery will influence in an undetermined 
fashion experimental capacities at low adverse tilts. As discussed at a 
later point, the effect of this pool on effective heat pipe length and ad- 
verse tilt were empirically determined in auxiliary tests. Heat transfer 
data was taken at adequate tilt to ensure recession of the liquid into the 
condenser. 

Experimental Technique 

Once the component parts of the heat pipe were cleaned and assembled, the 
heat pipe was flushed with several anhydrous ammonia charges and evacuated. 

The final charge used for data acquisition was approximately a 40% overcharge 
relative -to the arterial structure. This overfill was used to ensure priming 
of the axial artery. The influence of the pool resulting from overcharge on 
heat pipe effective length and adverse tilt was empirically determined in a 
glass tube of identical inside diameter as the heat pipe using the multitube 
artery and isopropyl alcohol to simulate ammonia priming behavior. After 
charging, the heat pipe- was allowed to remain stagnant for 1-2 hours in a hori- 
zontal position to aid in the venting of gas occlusions. At this point, the 
heat pipe was considered primed and charged. 



HEAT FLUX (WATTS) 
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EXTERNAL ADVERSE TILX CM 

FIGURE 3-5. PREDICTED CAPACITIES FOR THE THREE EVAPORATIVE HEAT TRANSFER 
SYSTEMS, INCLUDING VAPOR MICROCHANNEL PRESSURE DROPS 
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4.0 EXPERIMENTAL RESULTS 

The heat pipe evaporators were operated at a number of heat fluxes and eleva- 
tions to characterize response to heat flux density and liquid/vapor pressure 
differential. In this report, only the results from the 12.7/cm and 37.8/cm 
tests will be discussed in detail. Tests performed on the 25/cm surface are 
discussed in earlier reports. Figure 4.1 presents typical data of 
evaporator temperature difference versus heat flux for the 12.6 and 37.8 

groove/ cm surfaces at an adverse evaporator-condenser elevation of 6.4 cm. 

2 

At a heat flux density of 10 W/cm-, for example, the evaporator temperature 
drops are about 5.8° and 7.7°C, respectively, for 37.8 and 12.6 grooves/cm. 

Two distinct series of tests are shown for the 37. 8/ cm surface, and they 
demonstrate the importance of maintaining adequate groove-to-wick cover con- 
tact. The data shown with open faced triangular symbols was taken with an 
evaporator that had an estimated 0.002 to 0.004 cm diametrical clearance 
between the screen-covered mandrel and the grooved inner diameter. When this 
was changed to a diametrical interference fit of 0.002 to 0.004 cm, the data 
shown as solid-faced triangles were obtained. At high heat fluxes, the tight- 
fitting wick displayed significantly lower evaporator temperature differences, 
and it is assumed that the liquid film was prevented from receding into the 
screen, as was possible with the loose-fitting system. 

Figure 4.2 presents experimental data on heat transfer coefficients for the 

12.7/cm and 37.8/cm surfaces measured as a function of heat flux density 

under conditions ranging from reflux (-0.25 cm external tilt) to 6.4 cm 

2 9 

adverse elevation and from a heat flux density of 0.3 W/cm to 20 W/cm , In 
all cases, the heat transfer coefficient is based on the internal diameter 
and length of the screw thread zone and on the average evaporator temperature 
as determined by the 12 evaporator thermocouples. Table 4.1 presents previous 
experimental data for the 25/cm surface in tabular form. 

Although there were significant variations in behavior related to heat pipe 

tilt and groove density, several generalizations can be made. First, at high 

adverse tilt both surfaces displayed a cube root dependence of heat transfer 

coefficient on heat flux density. A least-squares power law fit to data 
2 2 

between 0.5 W/cm and 20 W/cm results in 
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TABLE 4-1. SUMMARY OF PREVIOUS HEAT TRANSFER 
MEASUREMENTS ON A 25/CM INVERTED 
MENISCUS SURFACE 


Heat Flux 
Density 
{W/cm2) 

Heat Transfer 
Coefficient 
(W/cm^K) 

Vapor Core 
Temperature 

(oc) 

External Tilt 
(cm) 

1.33505 

1.05122 

19.35 

0.0 

1.33513 

0.890088 

19.1 


2.64488 

1.3704 

20.2 


3.9878 

1.7414 

21.3 


5.92633 

1.8934 

22.7 


10.3716 

2.81074 

26 


1.33063 

0.65227 

19.7 

1.46 

2.63774 

0.912712 

21.05 


3.98101 

1.36805 

22.4 


5.92053 

1.85016 

24.1 . 


8.37155 

3.17104 

25.4 


10.3221 

3.05386 

27.2 


12.599 

2.83124 

29.3 


15.2481 

2.65184 

31.9 


20.0585 

2.66735 

36,5 


1.3317 

0.752371 

19.7 

3.0 

1.33114 

0.735439 

19.8 


2.63809 

0.995506 

21.2 


3.98026 

1.42152 

22.7 


5.91701 

1.. 60353 

24.5 


8.36851 

2.52825 

25.5 


8.36629 

2.02084 

26.9 


10.3191 

2.76651 

27.6 


10.353' 

2.03001 

29.3 


12.5968 

2.62434 

29.5 


15.2458 

2.57531 ■ 

32.3 


20.0585 

2.66735 

36.5 
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TABLE 4-1 

(Continued) 


Heat Flux 
Density 

(W/cm2) 

Heat Transfer 
Coefficient 
(W/cm2K) 

Vapor Core 
Temperature 
(OC) 

External Tilt 
(cm) 

1.33083 

0.74348 

19.9 

6.10 

2.63737 

1.00281 

21.4 


3.97853 

1. 30873 

22.9 


5.91283 

1.7386 

25.8 


8,36321 

1.93593 

27.5 


10,3488 

1.93075 

30.1 


1.33229 

0.938232 

19.9 

9.30 

2.63971 

1.29398 

21.4 


3.98042 

1.55485 

22.9 


1,33146 

0.870236 

20.0 

12.1 

2.63857 

1.18321 

21.5 


1.33051 

0.751702 

20.0 

15.1 
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FIGURE 4-1. REPRESENTATIVE INVERTED MENISCUS HEAT TRANSFER BEHAVIOR 
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FIGURE 4-2. MEASURED HEAT TRAIMSFER COEFFICIENTS FOR SCREEN COVERED GROOVES WITH AMMONIA WORKING FLUID 
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H = 0.593 q 

a 

and (W/cm^K) (4-1) 

H = 0.821 q 

d 

o 

where is the surface heat flux density in W/cm , and the equations are . 
respectively for 12.6/cm and 37. 8/ cm surface groove densities. 

Second, both surfaces were operable up to heat flux densities in excess of 

20 W/cm . At 20 W/cm , the film temperature’ drops were about 12.7°C (H = 

1.57 W/cm^K) and 8.7°C (H = 2.3 W/cm^K) for groove densities of 12.6/cm and 

37.8/cm. Third, both surfaces showed heat transfer coefficients from 1 to 
2 

2.5 W/cm K under refluxing conditions. Under adverse tilt conditions, the 
25/cm inverted meniscus surface shows the characteristic .increase in heat 

O 

transfer coefficient with heat flux density up to about 8 W/cm . At. this 
point, the heat transfer coefficient began to gradually decrease. The peak 
heat transfer coefficient ranged from 2 to a maximum of 3.2, with a higher 
heat transfer coefficient being observed at lower adverse tilts. The 25/ cm 
inverted meniscus surface displayed the highest heat transfer coefficient, 
but it may be more directly related to the semicircular groove shape than 
to an optimum groove density. 

A summary of heat transfer results is given in Table 4-2, along with cor- 
responding data for open-groove operation with the 25/cm surface. This table 
shows that the covering of .circumferential grooves with a capillary wick 
material increased the maximum heat flux density by a factor of 5 to 7 and 
the heat transfer coefficient by a factor of three to four for the same 
thread density. The 12.6/cm capillary wick-covered surface also outperformed 
the 25/cm uncovered grooves. At the 4 W/cm maximum heat flux attained with 
the open-faced 25/cm surface, the 12.6/cm inverted meniscus surface has - 

demonstrated a heat transfer coefficient of 0.9 to 1.0 W/cm K in comparison. 

2 

to only 0.44 to 0.79 W/cm K for the 25/cm open groove evaporator. 
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TABLE 4.2. SUMMARY OF HEAT TRANSFER RESULTS 


Parameter 


Value 



Groove density, cm~^ 

12.6 

25 

25 

37.8 

Surface type^^^ 

IV 

IV 

OG 

IV 

Groove shape 

65° V-shape 

.cup-shape 

cup-shape 66° 

V-gro( 

Heat transfer « 

coefficient, W/cm'^K^^' 

1.57 

■ 2 to 3 0 

.44 to 0.79 

2.35 

Maximum heat' flux 
density, W/cm^ 

•30' 

20 ■ 

4 

20 

Maximum prime height, cm 

— 

15 

6.8 

— 

Vapor, core temperature 
range, oc 

15-35 

20-40 

20-40 

15-35 


(T5 ; 

' = inverted meniscus, OG = open groove 

^^Vhese are values at 20 W/cm2 - refer to Figure 4-2 for variation with 
heat flux density. 
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5.0 DATA INTERPRETATION 

It has been found that wick-covered grooves display relatively high heat 
transfer coefficients and high surface heat flux capabilities compared to 
conventional screw-thread grooves, and that the heat transfer coefficient 
furthermore increases as the cube root of heat flux density. The overall 
heat transfer coefficient has also been found to vary slowly with thread 
density, so that it is possible to design evaporators with relatively coarse 
threads and yet maintain satisfactory heat transfer characteristics. A 
factor of 3 increase in groove density results in approximately a 40% in- 
crease in heat transfer coefficient . If the dependence of heat transfer 
coefficient on groove density is assumed to follow a power law, the propro- 
tionali'ty relationship is approximately 

H a (5-1) 

where N is the linear groove density. This behavior is at variance with 
established models of meniscus heat transfer, since the major fraction of 
heat is generally considered to pass through a relatively abbreviated zone 
at the fillet tip, hence a doubling of groove density should result in a 
corresponding factor of two increase in overall heat transfer coefficient. 

It has not been possible to definitively establish the heat transfer mecha- 
nisms responsible for the performance characteristics of the inverted meniscus 
surface, but several factors indicate that nucleate boiling may play an 
important role. 

In Figure 5.1, heat transfer .behavior for the inverted meniscus, is compared 

with the 25°C pool boiling correlation of Fritz^^^ and the -33°C boiling 

7 13 \ 

data of Grottzman, et-al.' The latter data is for the Union Carbide high 
flux porous boiling surface. The inverted meniscus data very nearly parallels 
the pool boiling curve and the data of Grottzman, et al. Both boiling 
relations follow approximately a 1/4-root dependence of ”H" on heat flux 
density. At a given heat flux density, the two wick-covered grooved surfaces 
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respectively show a constant 56% and 220% performance increase relative to 
pool boiling, but a reduction of about 80% relative. to the highly effective 
porous boiling surfaces. Together, the data of Fritz and Grottzman show that 
significant changes can be made in surface structure and operating pressure 
without influencing the functional dependence of boiling heat transfer coef- 
ficient on heat- flux density, and that the inverted meniscus surface has a 
functional dependence very similar to pool boiling. 

Liquid Film Boiling 

If boiling occurs during operation of the inverted meniscus surface, then 
certain criteria must be met. to assure that vapor bubbles will propagate from 
nucleation sites located under the liquid fillets. Since the spiral grooves 
in the existing surfaces were machine cut, it is probable that the grooved 
surfaces contain a large population of potential nucleation sites of varying 
radius. The problem is complicated to a degree because the liquid film is 
not constant in thickness and the groove wall beneath the film is not constant 
in temperature. At the fillet tip, the liquid film is very thin, and the 
wall temperature is depressed because of constriction resistance. Near the 
root of the fillet, the heat transfer will be a minimum and the groove wall 
temperature will approach the exterior evaporator wall temperature, T^g. The 
thick liquid film and high wall temperature at the meniscus root are both 
conducive to boiling initiation. 

To activate a nucleation site, the temperature field in the liquid film must 

(14) 

allow the vapor bubble to form a hemispherical shape. Bergles and Rohsenow'’ ' 
have developed a model for the prediction of incipient flow boiling that may 
be used here to estimate the critical nucleation radius for meniscus boiling. 
The physical situation is shown in Figure 5-2. If the temperature profile in 
the liquid film is linear and the film surface isin saturated equilibrium 
with the vapor phase, then it is postulated that boiling will be initiated 
when the local saturated vapor pressure at y = r^^. exactly compensates for 
the surface tension pressure 2a/r^^. . The radius r^^.| corresponds to the 
smallest cavity that can propagate vapor bubbles. Vapor bubbles at sites 
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With radii greater than r^^. will spontaneously grow, whereas bubbles at sites 
with a smaller radius will not grow in size. That is, larger bubbles do not 
encounter any unfavorable ratio of surface tension pressure to differential 
vapor pressure. This criteria does not in itself guarantee boiling will occur, 
since this depends on the size distribution of nucleation sites and the condi- 
tions within each site. Rather, boiling could occur for sites having a radius 
greater than r^^^. . It will be shown that this criteria for boiling applies to 
thin films if r^^ is small compared to the film thickness. 

For a small film temperature drop, the vapor pressure rise above saturated 
film conditions at y = r^^ can be obtained from the differential form of the 
Clausius-Clapeyron equation,. 

'sr-Pys “ -<5-2) 


In the absence of noncondensables, this vapor pressure rise must be greater 
than or equal to the pressure required to maintain the bubble radius. The 
pressure rise within the bubble must be (referring to Figures 2.4-1 and 5-1) 


P 


b 




(5-3) 


where r is the meniscus local radius of curvature. 

The ratio of expression 5-2 to 5-3 is the ratio of available differential 
vapor pressure to surface tension pressure. If T^j^^ is obtained from the 
assumed linear temperature profile, this ratio is 



■b 



(1) . '~b 


p+ 


2 - r,,/r 


(5-4) 



5-6 


where <f> = —r -9 — — (5-5) 

^fg " Tvs) 

It can be shown that the two pressure differences are equal for critical radii 
of 


m 


d 

2 



- 


1 8 <f./d 

1* J 

<1- , 


(1 + i)j^^-/r r 


(l+<t„,/r) (5-6) 


For most room temperature organic and inorganic working fluids the factor * 

is quite small, and the characteristic roots represent either a nucleation 

radius r„.« d or a value slightly less than d. For example, the factor 

at 25°C for ammonia, with a 1°C wall superheat, is approximately 6.6 x 10 

cm. This also means that the factor (f)/r can be neglected for a fillet radius 

-4 

greater than or equal to approximately 10 cm. A fillet radius of this 
magnitude may be present near the meniscus tip, but certainly not at the 
meniscus root for typical heat pipe wick materials. Therefore, the terms 
containing r may be set to zero. 


If a hemispherical vapor bubble of radius r^^ is generated at the liquid- 
metal interface, this bubble must grow in size and depart the surface for 
boiling to be truly initiated. If the hemispherical bubble cap is considered 
to be a "Type One" geometry, and the larger bubbles shown in Figure 5.2 to 
be "Type Two" geometries, then the bubble will definitely grow if the pres- 
sure ratio P"’" at the bubble cap is always greater than or equal to 1.0. 
Neglecting the influence of r , the pressure ratio ?■*" for the Type Two 
geometry is 


p+ 


1 - r,^/d (1+a)’ 




1 

2<l> 


(5-7a) 
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and where r^^^. is given by the smaller root of equation 5-6, with set 

to zero. If the ratio p+ is calculated as a function of r^ and r^^^. , it is 
found that for >"p^/d ^ 0.17, the Type One hemispherical cap bubble will 
spontaneously grow into a Type Two bubble. In contrast, for r^^ s 0.17, the 
final bubble height equals the initial radius; that is, y^ = r^^. . Figure 5.3 
shows the final size of the bubble, y^/d, as a function of ^'p,^-/d. 

It is worthwhile to note that with this model for incipient nucleate boiling, 
the vapor bubble never actually penetrates the surface. However, it can be 
assumed that a 90% to 95% penetration is adequate to breach the film, since 
it can be expected that the actual vapor pressure within the bubble is some- 
what higher than the characteristic pressure at its extreme tip. Therefore, 
for »'p^/d ^ 0.05, Figure 5.3 shows that the critical state for boiling 
initiation is the production of a Type 1 bubble cap, and the relation between 
nucleation site radius and wall superheat at this critical condition is 

(boiling Initiation (5-8) 

when r„. ^ 0.05 d) 
m 

For nucleation site radii greater than 0.05 d, additional superheat is needed 
to attain a 90% to 95% film penetration. If the superheat factor necessary 
to form a bubble of height y = Bd is denoted then the ratio of to 
^yf s, for constant film thickness. 



where r+^. = 




0.05 


0.10 


0.15 


O.ZO 


F5GURE 5-3. THE EQUILIBRIUM SIZE FOR AN INITIALLY HEMISPHERICAL VAPOR BUBBLE OF 
RADIUS r„j, SHOWING THE GROWTH PHENOMENA CHARACTERISTIC OF NU- 
CLEATION AT SMALL SITES 
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This ratio is presented in Figure 5-4 for B = 0.50, 0.75, 0.90, and 0.95. 

As an example, for a bubble that initially in Type 1 condition spans 20% 
of the film thickness = 0.20), growth of the bubble requires additional 
wall superheat as follows. 


For y^ = 0.5 d; 

= 0.75 d; 
= 0.90 d; 
= 0.95 d; 


'^nfV 


1.1 

1.6 


3.4 

6.4 


At a constant film thickness, a hemispherical bubble that is initially 20% 

of the thickness will require the wall-to-vapor temperature difference to 

increase by a factor of 3.4 to enlarge the bubble to 90% of the film depth. 

In the particular case of ammonia at 25°C, let r . = 2.5 {10"^) cm, and 
—3 

d = 1.25 (10 ) cm. The superheat corresponding to the initial condition is 

0,33°C, while for y^ = 0.9 d, the superheat required is 1.1°C. Although this 
analysis is necessarily approximate, it indicates that vapor bubbles can be 
produced within the liquid fillets in a grooved or inverted meniscus heat pipe 
at relatively low wall superheats, and that the bubbles may be quasi -stable 
over certain heat flux ranges. These stable bubbles could lead to enhanced 
heat transfer through capillary jet phenomena. 

For -the "H" versus q data gathered in this experimental program, the incipi- 
ent nucleation radii can be calculated at each condition using the equation 


(r+. < 0.05) (5-10) 

To establish for the present work the nucleation site radii that might be 

activitated, and the heat fluxes necessary, the characteristic film thickness 

d appearing in Equation 5-10, has been set equal to the dimension d of 

-3 ^ 

Figure 5-2. This value is approximately 1.32 x-'lO cm. Predicted values of 
2r^^ are presented in Figure 5-5. 


r: 


m 



« ° ''fc Tys H 
^ • '’fg \ 





DJMENSIONLESS CAVITY SIZE. r„|/d 
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FIGURE 5-4. THE ADDITIONAL SUPERHEAT REQUIRED TO FORM A VAPOR BUBBLE OF HEIGHT Yp FROM A HEMISPHERICAL BUBBLE 
CAP OF RADIUS 
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FIGURE 5^. THE MINIMUM PORE RADIUS (r„j) FOR INCIPIENT BOILING, AND THE COR- 
RESPONDING VAPOR BUBBLE DEPARTURE RADIUS, 
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For most machined surfaces, nucleation sites range in size from 10” cm to 10” 

cm-. These calculations indicate that all incipient nucleation site diameters 

-5 -4 

fall between 3 x 10 cm and 3 x 10 cm, which is well within the range of 
observed nucleation site sizes. 


Vapor Bubble Departure Radius 

For spacecraft applications- it is important to determine if nucleate boiling 
characteristics' in a thin, surface tension dominated liquid meniscus will 
change in going from Ig to zero g conditions. One of the most fundamental 
differences, would be the possible formation of discrete vapor bubbles in Ig 
only, as a result of interacting surface tension and bouyancy forces. The 
departure radius for a vapor bubble is given by the approximate transcendental 
formula^^®^ 


, 3^3 
bi 2 


°-^ni 



"bi^'ni 


(5-11) 


Figure 5-5 shows the critical nucleation diameter and corresponding bubble 
departure diameter for the "H" versus q data of Figure 4-2. The bubble 
departure diameters fall in the range 1.5 (10” ) to 4 (10” } cm, and are much 
larger than the characteristic film thickness, d^, which is 1.32 (10”^) cm. 
Since the lift-off diameters based on bouyancy exceed the film thickness, it 
can be expected that gravity will play a minor role in any meniscus boiling 
phenomena, and that zero-g and 1-g performance should be very similar. 


Original page is 
OF POOR QUALITY 
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6.0 SUMMARY 


Inverted meniscus surfaces of 12,6/cm and 37. 8/ cm V-groove density have 

2 

been operated to heat flux densities in excess of 20 W/cra with ammonia. ' 

2 2 
At 20 W/cm , the respective heat transfer coefficients were 1.57 W/cm K 

2 

and 2.3 W/cm K, while at lower heat flux densities, the overall evaporative 
heat transfer coefficient varied as the cube root of heat flux. If the 
heat transfer coefficient is assumed to follow a power law dependence on 
groove density, then the power is +0.32. That is, the heat transfer coef- 
ficient changes very slowly with variations in groove density, and the 
inverted meniscus is of particular interest .in that relatively high heat 
transfer coefficients and very high heat flux densities have been demonstrated 
for low groove densities. 

The approximate cube root dependencies on heat flux and groove density are 
unlike the dependencies normally expected for grooved evaporative surfaces, 
where the heat transfer coefficient is assumed to be linear in groove density 
and relatively independent of heat flux. Instead, the surfaces show behavior 
similar to pool boiling. 

Pool boiling heat transfer coefficients for ammonia follow an approximate 
fourth-root dependence on heat flux density under diverse surface conditions. 
Furthermore, calculations show that conditions are favorable for initiation 
of thin-film boiling within the liquid fillets for a wide range of heat flux 
densities. These calculations further imply that nucleate boiling may be 
present in open groove evaporators as well as in inverted meniscus evaporators 
and that the bubbles produced may influence both heat transfer and fluid trans 
port. 

To gain a fundamental understanding of fillet evaporation, heat transfer 
models for grooved surfaces were developed. Two-dimensional conduction models 
for rectangular grooves showed that thermal constriction resistance in each 
metal tooth is a negligible factor in restricting evaporation of ammonia for 
rectangular grooved metal surfaces constructed of high conductivity metals 
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such as aluminum. With low conductivity metals such as stainless steel, 
however, this resistance could be quite significant. 

A composite model for evaporation from rectangular grooves was also developed 
that includes one-dimensional conduction heat transfer in the metal tooth 
and liquid fillet, and an approximate solution to the Navier Stokes- equation 
for liquid flow within the fillet. In the past, investigators have- modeled 
meniscus heat transfer as a pure two-dimensional conduction problem. The 
present analysis indicates that liquid flow from the base of a fillet to 
the tip can create a significant liquid pressure drop parallel to the groove 
wall, and that this liquid pressure change can strongly influence fillet 
shape and thermal resistance. These calculations suggest that the modeling 
and optimization of grooved surfaces for evaporative heat transfer may 
require consideration not only of conduction processes, but also of fluid 
transport processes within the liquid menisci. 
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